Scaling Behavior of Entanglement in Two- and Three-Dimensional Free Fermions 
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Exactly solving a spinless fermionic system in two and three dimensions, we investigate the scaling 
behavior of the block entropy in critical and non-critical phases. The scaling of the block entropy 
crucially depends on the nature of the excitation spectrum of the system and on the topology of 
the Fermi surface. Noticeably, in the critical phases the scaling violates the area law and acquires a 
logarithmic correction only when a well defined Fermi surface exists in the system. When the area 
law is violated, we accurately verify a conjecture for the prefactor of the logarithmic correction, 
proposed by D. Gioev and I. Klich ^quant-ph/050415]j. 
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The nature of many-body entanglement in various 
solid-state models has been the focus of recent interest. 
The motivation for this effort is two-fold. On the one side, 
these systems are of interest for the purpose of quantum 
information processing and quantum computation 0] . At 
a fundamental level, the study of entanglement represents 
a purely quantum way of understanding and characteriz- 
ing quantum phases and quantum phase transitions in 
many-body physics 

A striking feature of entangled states \'^) is that a 
local accurate description of such states is impossible, 
namely each subsystem A of the total system U can 
have a finite entropy, quantified as the von- Neumann en- 
tropy Sa — — TrpA log2 PA of its reduced density matrix 
PA — Tr[/\^ I 1, whereas the total system clearly has 
zero entropy. The entropy of entanglement Sa of the 
subsystem is a reliable estimate of the entanglement be- 
tween the subsystem A and the rest, U \ A. Assuming 
that the system U corresponds to the whole universe in d 
dimensions, a fundamental question concerns the scaling 
behavior of the entropy of entanglement Sl of an hyper- 
cubic subsystem L'^ (hereafter denoted as a block) with 
its size L. Indeed, such scaling probes directly the spatial 
range of entanglement: when the block size exceeds the 
characteristic length over which two sites are entangled, 
the block entropy should become subadditive, and scale 
at most as the area of the block boundaries, following a 
so-called area law: Sl ~ L'^~^. A crucial question is then 
if and how the scaling of the block entropy changes when 
the nature of the quantum many-body state evolves in a 
critical way by passing through a quantum phase transi- 
tion, and how the characteristic spatial extent of entan- 
glement relates to the correlation length of the system. 

This question has been extensively addressed in the 
case of one-dimensional spin systems [a,IESIl|j in chains 
of harmonic oscillators | 9| , llOl and in related conformal 
field theories (CFT) Here it is found unam- 

biguosly that in states with exponentially decaying (con- 
nected) correlators Sl follows the area law Sl ^ LP ^ 
saturating to a finite value, whereas for critical states, 



displaying power-law decaying correlations, a logarith- 
mic correction to the area is always present: Sl = 
[(c -I- c)/6]log2L, where c is the central charge of the 
related CFT. The asymptotic value of the block entropy 
is found to diverge logarithmically with the correlation 
length, 5*00 ~ log2(C), which clearly establishes the re- 
lationship between entanglement and correlations. The 
above picture holds true only in presence of short-range 
interactions; on the contrary, in presence of long-range 
interactions the divergence of the correlation length can 
be still accompanied by the area law 0, 0| . 

In higher spatial dimensions less results are available, 
and the general relationship between the block entropy 
scaling and the correlation properties of the quantum 
state is still unclear even for short-range interactions. 
In free-boson systems, it has been generally proven that 
the area law is satisfied for non-critical systems 0, ■ 
For free-fermion systems, on the other hand, it has been 
proven U M that critical systems with short-range 
hoppings and a finite Fermi surface exhibit a logarith- 
mic correction to the area law 

SL = C/?,{\og^L) L^-K (1) 

In this paper, we investigate a general quadratic 
fermionic Hamiltonian both in d = 2 and d = 3. Upon 
tuning the Hamiltonian parameters, this model has dis- 
tinct critical phases with and without a finite Fermi sur- 
face, as well as non-critical phases. The scaling behavior 
of the block entropy is accurately obtained through ex- 
act diagonalization. In non-critical states we find that 
the area law indeed holds, and we confirm that logarith- 
mic corrections to such law are present in critical states 
with a finite Fermi surface, as found in Refs. 13 [l^ The 
prefactor C of the L-dependence of Sl in Eq. is found 
to be very accurately pred icted by a formula based on 
the Widom conjecture [l^. On the other hand, for crit- 
ical states with a Fermi surface of zero measure, we find 
that the corrections to the area law are either absent or 
sublogarithmic. This means that the relationship between 
entanglement and correlations in higher dimensional sys- 
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FIG. 1; Phase diagram of the model Eq.Q for the case d — 
2. The roman numbers for the various phases are explained 
in the text. Representative contour plots of the dispersion 
relation Afc are also shown. There the black areas corresponds 
to Afe = and the white areas to the top of the band. 

tems is different than in d = 1, and that a crucial role 
is played by the geometry of the Fermi surface or, al- 
ternatively, by the density of states at the ground state 
energy. 

We consider a bilinear spinless fermionic system on a d- 
dimensional hypercubic lattice with hopping and pairing 
between nearest-neighbor lattice sites 

H = J2 Re.- " 7(44 + ^^•^^)] - E 2A4ci. (2) 

A is the chemical potential, while 7 is the pairing po- 
tential. The sum of extends over nearest- neighbor 
pairs. The above Hamiltonian is a d > 1 generaliza- 
tion of the Id spinless fermionic Hamiltonian which is 
obtained by Jordan- Wigner transformation of the XY 
model in a transverse field Although in d > 1 
the exact relationship to the XY model is lost, we can 
imagine the above Hamiltonian to represent the effec- 
tive fermionic degrees of freedom of an interacting system 
with quantum-critical phases. 

A more insightful expression for the Hamiltonian of 
Eq. (0) is obtained upon Fourier transformation to mo- 
mentum space: 

H = [^-2tfcC^Cfc -I- iAfc(4c^fc + c-kCk) 

k 

d d 

tk = — COS ka Afc = 7 sin ka (3) 

a— 1 a— 1 

The pairing potential in fc-space, Afc, clearly reveals a 
p-wave symmetry. 



This Hamiltonian can be diagonalized exactly by a Bo- 
goliubov transformation, to give 

H = Y.Kkfih Ak = 2^tl + Al (4) 

k 

Depending on the parameters 7 and A, this system has a 
rich phase diagram, including metallic, insulating and (p- 
wave) superconducting regimes, as shown in Fig. ^ The 
different phases are certainly distinguished by the differ- 
ent decay of the correlation function, which tells apart 
the critical from the non-critical phases. Nonetheless, a 
classification which turns out to be relevant for the study 
of entanglement is based on the features of the gapless 
excitation manifold Afc = 0. Such manifold can be char- 
acterized by the density of states at the ground-state en- 
ergy (7(0), and by the so-called co-dimension |l8[ Il9j| d, 
defined as the dimension of fc-space minus the dimension 
of the Afc = manifold. We notice that the existence of 
a finite Fermi surface at zero energy implies that d = 1 
and g(0) > 0, while in absence of a finite Fermi surface 
we have d > 2 and g{0) > or 17(0) = depending on 
the dispersion relation Afc around its nodes. 

According to this classification, which turns out to be 
relevant for the study of entanglement, we can distinguish 
three phases: 

• Phase /, {0 < A < d, 7 = 0}, and {A = 0, 7 > 0} if 
d — 2. For 7 = 0, Eq. jSJ reduces to a simple tight- 
binding model, which is in a metallic state with a 
2d-io\d symmetric Fermi surface as far as A < d. 
In d = 2, for A = the system is still a metal 
with a well defined Fermi surface, which is simply 
kx ~ fcy ± TT, and whose symmetry is lowered by 
the presence of the 7 term in the Hamiltonian. In 
this phase, g{0) > 0, and d = 1 everywhere except 
at the point {A = d, 7 = 0} where d = 2. 

• Phase //, {0 < A < d,7 > 0}, and {A = 0,7 > 0} 
if d = 3. Away from the boundary lines of this 
phase, the system is in a p-wave superconducting 
state, with a finite pairing amplitude (4cl^) 7^ 0. 
Such pairing amplitude vanishes at the boundaries 
of this region. The dispersion relation Afc has point 
nodes in d = 2 and line nodes in d = 3. Everywhere 
in this phase g{0) ~ and d ~ 2. 

• Phase III, {A > d}. In this phase the system is 
in an insulating state with a gap in the excitation 
spectrum. Here g(0) = and d = d. 

This shows that, in terms of the spectral properties, 
the above system has two distinct critical phases, {I and 
//), which are both gapless, and a non- critical phase 
{III}- Numerical evaluation of the correlation function 
through integration over the first Brillouin zone (FBZ), 

yFBz(27r)'^2Afc' ' 
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shows an expected power-law decay in the critical phases 
and an exponential decay in the non-critical one. 

We then proceed to the evaluation of the block entropy 
of entanglement. The ground state of Eq.® is known to 
be a Gaussian state, whose density matrix can be ex- 
pressed as the exponential of a quadratic fermion oper- 
ator [23, 01 . To obtain the reduced density matrix of a 
L''- subsystem, Grassman algebra is needed . Using a 
Bogoliubov transformation, the reduced density matrix 
PL can then be written as 



where df ,di are the new Fermi operators after the trans- 
formation, and ^ is a normalization constant to en- 
sure Tr{p) = 1. The single-particle eigenvalues e; can 
be obtained from {c'^Cj) and {cfc'^) by the following 
formulaji^ ■ 



where dj = (clcj) and Fij — {clc1j)]P is the orthog- 
onal matrix that diagonalizes the left side of the above 
equation. The Block entropy can then be calculated in 
terms of ei as: 



In d = 1 the above formulas reproduce the scaling 
of the block entropy as observed in the XY model in 
a transverse field [^0. In d — 2 the phase diagram is 
richer, and we need to consider the various phases one 
by one. We begin with the critical metallic phase (I), 
7 = 0, < A < d. For this case a logarithmic correction 
to the area law, Sl — (C(A)/3)(log2 is observed 

for all values of A, as shown in Fig. El This is in full 
agreement with the results of Refs.|l5L Il6l|. which pre- 
dict this behavior in presence of a finite Fermi surface. 
More specifically, Ref. also supplies us with an ex- 
plicit prediction for the A dependence of C(A), based on 
the Widom's conjecture jl^, in the form 

C{X) = / / ■ np\dS^dSp (9) 

where is the volume of the block normalized to one, 
r(A) is the volume enclosed by the Fermi surface, and the 
integration is carried over the surface of both domains. 
A numerical fit of the calculated asymptotic behavior of 
Sl through the formula Sl = ^L'^~^ \og^{L) + BL"^-^ + 
^ provides us with the exact result for the C(A) 
prefactor. In Fig. |3|the prediction of Ref. 0, Eq.®, 
for the case {0 < A < d, 7 = 0} is compared to our 
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FIG. 3: A-dependence of the C coefficient in Eq.Q in d = 2 
and d = 3. The values extracted from fits to our numerical 
data are compared with the predictions of Ref. IT^ In d = 2, 
the exact form of C(A) can be obtained, which is equal to 
lcos-\X-l) 



numerical results both for d — 2 and d = 3. The agree- 
ment is clearly striking. Moreover, for {A = 0,7 > 0} in 
d = 2 the formula Eq.@ predicts C — I, which is also 
accurately verified by our data as shown in Fig. El This 
proves that the formula Eq.© is essentially providing 
a complete analytic form for the leading behavior of the 
block-entropy scaling in arbitrary dimensions for systems 



Pl = A exp - ^ejd/ 



(6) 



(^C - F - ^){C + F - ^) = FIG. 2: Scahng of the block entropy & in d = 2 for 7 = (left 

2/x 2 / \ 2/'M panel) and A = (right panel) . The solid lines correspond to 

Pdiagjtanh (^),tanh (f),...,tanh (^)|p-i(7) fits according to the formula 5i = fLlogjCL) A. 



--?{ln[l + exp(-.0] + ^5^^} (8) 
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Further investigations in systems with d > 1 are clea rly 
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FIG. 4: Scaling of the block entropy Sl in d = 2 for 7 
(left panel) and A = 1 (right panel). 



with a finite Fermi surface. 

We then turn to the other two phases, // and ///. Two 
scans through these phases, at fixed 7 = 1 and at fixed 
A = 1 are shown in Fig. 0] We observe that logarithmic 
corrections are absent in both, and only sublogarithmic 
corrections are possible. For {A > d, 7 = 0} 5l = iden- 
tically, and the state is not entangled. While the area law 
is expected to hold in the non-critical phase III, it is sur- 
prising to observe it enforced also in the critical phase //, 
which has a diverging correlation length. This clearly re- 
veals that the connection between block-entropy scaling 
and correlation properties is not as straightforward as in 
d > 1. 

Our results for the entanglement behavior, co- 
dimension, density of states and correlation properties 
are summarized in Tabled A crucial difference between 
the two critical phases / and // is the co-dimension d, 
and the density of states at the ground state energy. We 
have d — 1 and g{0) > in the phase /, which shows 
logarithmic corrections to the area law, whereas d = 2 
and g{0) = in the phase //, in which the area law 
is verified up to sublogarithmic corrections. It is there- 
fore tempting to conjecture that a codimension 0? < 1 
or, alternatively, a finite density of states at the ground 
state energy g{0) > is a necessary condition for critical 
phases in d > 1 to show violations of the area law. For 
the fermionic system under consideration, d = 1 requires 
the existence of a finite Fermi surface, which is the ba- 
sic assumption of the proof of area-law violation in Refs. 
Iisl nil This conjecture would generalize the results for 
d = 1, where the co-dimension can only take the value 
d = 1, and only critical phases with g{0) > have been 
explored in the literature. 
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power-law decay 
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power-law decay 
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exp. decay 



TABLE I: Summary of the entanglement scaling properties, 
co-dimension, density of states and decay of correlations in 
the three phases of the model Eq.Q in d= 2, 3. 



needed to confirm this picture, and to clarify whether 
more severe violations of the area law are possible in pres- 
ence of infinitely degenerate ground states or in systems 
with a fractal Fermi surface .15i|. During the c omp letion 
of this manuscript we became aware of Ref. whose 
results are in full agreement with the ones reported in 
our work. 
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